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Abstract 

The results of a numerical study of the magnetic dynamo effect in cylindrical von Karman 
plasma flow are presented with parameters relevant to the Madison Plasma Couette Experiment 
(MPCX). This experiment is designed to investigate a broad class of phenomena in flowing plasmas. 
In a plasma, the magnetic Prandtl number Pm can be of order unity (i.e., the fluid Reynolds 
number Re is comparable to the magnetic Reynolds number Rm). This is in contrast to liquid 
metal experiments, where Pm is small (so, Re » Rm) and the flows are always turbulent. We 
explore dynamo action through simulations using the extended magnetohydro dynamic (MHD) 
NIMROD code for an isothermal and compressible plasma model. We also study two-fluid effects 
in simulations by including the Hall term in Ohm's law. We find that the counter-rotating von 
Karman flow results in sustained dynamo action and the self-generation of magnetic field when the 
magnetic Reynolds number exceeds a critical value. For the plasma parameters of the experiment 
this field saturates at an amplitude corresponding to a new stable equilibrium (a laminar dynamo). 
We show that compressibility in the plasma results in an increase of the critical magnetic Reynolds 
number, while inclusion of the Hall term in Ohm's law changes the amplitude of the saturated 
dynamo field but not the critical value for the onset of dynamo action. 



I. INTRODUCTION 



The dynamo phenomena, where magnetic fields are self-generated by a moving and elec- 
trically conducting fluid, is one of the most intriguing subjects of modern magnetohydro- 
dynamics (MHD). Dynamos have particularly important applications in astrophysics [1, 2]. 
Today it is widely believed that the magnetic fields of planets and stars originate from 
dynamo action in their interiors [3-5]. As far as we know, all astrophysical dynamos are 
turbulent with extremely high fluid Reynolds numbers i?e, which makes theoretical treat- 
ments very involved and realistic simulations are currently impossible. At the same time, 
some key physical processes likely at work in astrophysical dynamos can be revealed by 
considering idealized laminar dynamos related to spatially smooth stationary velocity fields 
at comparatively low Re. In particular, the theoretical study of laminar dynamos allows one 
to determine the critical magnetic Reynolds number, above which dynamo excitation takes 
place, to find the structure and magnitude of saturated dynamo field, and to understand 
the influence on the dynamo of different plasma effects, such as compressibility, two-fluid 
effects, and anisotropic transport. 

A number of laminar MHD flows appropriate for dynamo generation have been analyzed 
in the literature [6-11]. In most of these studies only the kinematic dynamo problem has 
been considered, in which the magnetic induction equation is solved as an eigenvalue problem 
for a given velocity field (not necessarily satisfying the Navier-Stokes equation) to find the 
growth rate of the magnetic field. The non-linear feedback of the fields on the flows is ignored 
in the kinematic treatment. Among the laminar flows that lead to kinematic dynamos 
at sufficiently high magnetic Reynolds numbers are cylindrical helical jets (Ponomarenko 
dynamo [6]), two-dimensional spatially periodic arrays of helical jets (Roberts' scheme [7]), 
cylindrical von Karman [8] and Taylor-Couette flows [9], and spherical Dudley- James flows 
[10]. The first three of these flows have been tested recently in experiments with liquid 
sodium, and successful observations of dynamo action have been reported in Refs. [12- 
14]. The flows in these experiments were turbulent, thus making it difficult to compare 
experimental data with predictions of laminar dynamo theory, although in the first two, 
the role of turbulence was small as the flows were strongly constrained. This is a common 
disadvantage of all liquid metal dynamo experiments: the extremely low magnetic Prandtl 
numbers (the ratio of kinetic viscosity v to resistivity rj or, equivalently, the ratio of magnetic 
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Reynolds to fluid Reynolds Pm = v /rj = Rm/Re ~ 1CT 5 for liquid sodium) requires very 
high fluid Reynolds numbers (Re ~ 10 6 — 10 7 ) in order to achieve the magnetic Reynolds 
number sufficient for dynamo action (Rm ~ 10 1 — 10 2 ) in liquid metals. As a result the 
relevant flow is always turbulent. 

The present paper is motivated by the Madison Plasma Couette Experiment (MPCX) 
[15], which is designed to study MHD phenomena driven by plasma flows. One of the 
novelties of this experiment is the ability to change the magnetic Prandtl number of the 
plasma by several orders of magnitude from Pm « 1 to Pm » 1. This flexibility 
makes it possible to investigate laminar dynamos by choosing a regime with Pm ~ 1 — 10 
and Re ~ 10 2 . As a result the direct comparison of experimental data with numerical 
simulations of laminar dynamo can be performed. Such a comparison can also be used for 
the first time to test different MHD models as well as the numerical codes which simulate 
them. 

The goal of this study is to investigate possible dynamo action in the MPCX using the 
extended MHD code NIMROD [16], which can accurately model plasma dynamics in the 
specific geometry for realistic experimental conditions. Among the features of NIMROD is 
the possibility to study effects beyond standard MHD, including the addition of the Hall term 
in Ohm's law. The effect of the Hall term on dynamo action has recently been studied in 
periodic box simulations [18, 19], and the Hall term will almost certainly play an important 
role in MPCX. The results of our simulations can also be used for the optimization of plasma 
parameters and as a guidance for the experimental operation. 

In this paper we report the results of NIMROD simulations of laminar magnetic dynamo 
in the cylindrical von Karman flow under conditions relevant to MPCX. Simulations are done 
for an isothermal compressible MHD plasma model with and without two-fluid effects (the 
Hall term). The structure of the paper is following. In section II we briefly review the MPCX 
experiment and describe the NIMROD plasma model. In section III the hydrodynamical 
properties of von Karman flow are studied. In section IV the kinematic dynamo problem is 
considered for von Karman flow in a cylinder, and the self-generation of the magnetic field 
is demonstrated for parameters achievable in the experiment. In section V the results of 
simulations of nonlinear dynamo saturation are presented, and the effect of the Hall term 
is studied. In section VI we conclude with a discussion of the effects that can influence the 
magnetic dynamo in such plasma flows. 
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II. NIMROD MODELS OF MPCX 



MPCX is closely related to the spherical plasma experiment described in Ref. [20], though 
here the geometry is cylindrical and the apparatus is somewhat smaller in size (1 m in diam- 
eter). In MPCX the plasma is confined by a multi-cusp magnetic field created by axisym- 
metric rings of permanent magnets of alternating polarity and localized at the boundary of 
the cylindrical chamber (Fig. 1). There are 10 magnetic rings at the cylindrical wall and 8 
at both the top and bottom end-caps. Ring anodes and cathodes positioned between the 
magnet rings can be biased with arbitrary potentials. The resulting E x B drift of plasma is 
in the azimuthal direction and can be an arbitrary axisymmetric function at the boundary 
of the vessel. This arrangement allows arbitrary shear flows to be imposed in the MPCX 
experiment. 

The use of a plasma gives experimentalists great flexibility in choosing the regimes of 
operation. By varying the plasma density (gas flow rate), ion mass (H, He, Ne, Ar), electron 
temperature (heating power), and flow speed (bias potentials of electrodes), a wide range of 
parameters can be achieved in experiment (Table I). Such flexibility is advantageous over 
the liquid metal dynamo experiments, where Pm is fixed and small. For description of the 
plasma parameters in our simulations we introduce several standard dimensionless numbers: 

Magnetic Prandtl: 

V n [10 18 m- 3 ]/4 /2 A 2 

Fluid Reynolds: 

Re = Ml = o 52 ^HV [km/s]rio[10 18 m- 3 ]^ /2 A 
v ' 7f /2 [eV] 



Magnetic Reynolds: 



Rm = RoVo - 21 RoH V ° [km/S] Tfi3/2[eV] , (3) 
'/ A 



Mach: 

M ^ . 10 M< (4) 

C s 7 V2 Te 1/2 [eV] 

Hall: 

1/2 

£ = c = o 23 — (5) 

R^ pt ' i? [ m ]n; /2 [10 18 m- 3 ]' 
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permanent magnets 




FIG. 1: Madison Plasma Couette Experiment (MPCX): (a) sketch; (b) partial vertical cross section. 
Rings of permanent magnets of alternating polarity line the inside of the cylinder with their poles 
oriented normally to the walls. Ring anodes and cathodes are placed between the magnets. The 
resulting E x B drift is in the azimuthal direction. By varying the potential between the anodes 
and cathodes the velocity forcing at the outer boundary can be customized. 
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where v is the plasma kinematic viscosity, 77 is the magnetic diffusivity, A is the Coulomb 
logarithm (typically A ps 10 — 20), C s is the ion sound speed, 7 is the adiabatic index, c is 
the speed of light, u pi is the ion plasma frequency and other parameters are defined in Table 
I. Formulas (l)-(3) for numerical estimates of Pm, Re and Km are derived from Braginskii 
equations for a plasma with singly charged ions in a weak magnetic field [21]; the weak-field 
approximation is reasonable for the bulk of MPCX because the high-multipole cusp field 
is concentrated mostly near the wall and quickly falls off away from it. The typical values 
of these non-dimensional numbers are listed in Table I. For convenience we also give the 
"inverse" mapping formulary: 



Peak driving velocity, km/s: 

Pi K l m 



^0 = 2-54 ' ; , (6) 



Average number density, 10 m 



18 ™-3. 



Electron temperature, eV: 



Ion temperature, eV: 



n = 0.053 ^ (7) 
e 2 Rq [m] 



Te = 0.065 A ^ K ™ , (8) 
7 V4MV2 J R 1 / 2 [ m ]' 



T^0.35 Al/2 ^ 3/20 ^ 10M3/1 °. (9) 

£4/5^2/5^/2 [ m ] 



Results presented in this paper are obtained using the extended MHD code NIMROD 
[16]. As a simulation framework, we choose the isothermal Hall MHD approach. This is one 
of the simplest NIMROD models allowing for two-fluid effects and compressibility, and this 
appears to be a good approximation for the plasma under experimental conditions. The 
equations of this model in non-dimensional form are: 

^ = -V-(nv), (10) 

n ^ = " n ( v • V ) v - Jp Vn + (V x b) x b + — ( V 2 v + -V(V • v) J , (11) 

^ = Vxfvxb--(Vxb)xbV^V 2 b. (12) 
or V n J Km 
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TABLE I: Parameters of MPCX 



Quantity 


Symbol 


Value 


Unit 


Radius of cylinder 


Rq 


0.5 


m 


Height of cylinder 


H 


1 


m 


Peak driving velocity 


V 


0-20 


km/s 


Average number density 


no 


10 17 - 10 19 


m -3 


Electron temperature 


T e 


2-10 


eV 


Ion temperature 


Ti 


0.5-4 


eV 


Ion species 




H, He, Ne, Ar 




Ion mass 


Hi 


1, 4, 20, 40 


amu 


Ion charge 


z 


1 


e 


Magnetic Prandtl 


Pm 


1 x 10~ 3 - 5 x 10 2 




Fluid Reynolds 


Re 


- 4 x 10 4 




Magnetic Reynolds 


Rm 


- 1 x 10 3 




Mach 


M 


0-4 




Hall 


e 


0.15-1.8 





In these equations r, n, v and b stand for normalized time, number density, velocity and 
magnetic field, respectively: 

V 0+ p V B 

T — —t, n= , v = — , b= , 

Ro n mi V V y/A7rnomi 

where p is the mass density, rrii is the ion mass. The unit of length is the cylinder radius 
Rq, while Vq is the peak velocity of the driven von Karman flow. An important difference of 
this system of equations (10)-(12) from a standard single-fluid MHD model is the inclusion 
of the Hall term in the magnetic induction equation (12), which takes into account two-fluid 
effects. The magnitude of this term is characterized by the non-dimensional Hall number e. 
We will consider simulations where the Hall term is significant (0 < e < 1 ) and others where 
it is absent [e — 0). The simulations are performed in a non- rotating cylindrical coordinate 
system (r, ip, z), with the plasma occupying the region (0 < r < 1, — 1 < z < 1). 

The Hall MHD equations (10)-(12) also require the specification of boundary conditions. 
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Two different sets of boundary conditions are used in the simulations. Set I is used only to 
demonstrate the possibility of stirring the plasma with the applied multi-cusp magnetic field 
and an appropriately modulated tangential electric field at the boundary. In set I no-slip, 
stationary, rigid walls are assumed, so all components of the velocity vanish at the boundary, 

v| r = 0. (13) 

For the magnetic and electric fields in set I, we assume perfectly conducting walls, imply- 
ing that the time-varying normal component of the magnetic field and the time-varying 
tangential component of the electric field are zero at the boundary: 

b n \ T = 0, Et\r = 0, (14) 

where the normalized electric field is 

E = -v x b + -(V x b) x b + V x b. 

n Rm 

Note that the externally applied time-independent multi-cusp magnetic field b and bound- 
ary electric field E do not satisfy conditions (14), and thus provide E x b stirring at 
the boundary. Using different modulations of the tangential boundary electric field we have 
successfully simulated several types of flows. Fig. 2 shows the results for the so-called von 
Karman flow, in which plasma is driven in opposite azimuthal directions near the top and 
bottom end-caps. Having demonstrated that Eo x bo stirring will successfully drive a von 
Karman flow under realistic experimental conditions, we now turn to a simpler set of bound- 
ary conditions. 

All results reported in the later sections are obtained with the boundary conditions of 
set II. In set II we ignore the applied multi-cusp magnetic field and the tangential electric 
field. Instead we assume that the driving of the plasma is due to differentially rotating walls. 
This assumption greatly simplifies the model and allows us to focus on the physics of the 
dynamo action and not on details of the plasma driving. The boundary conditions for the 
full electric and magnetic fields are: 

6 n |r = 0, E t \ r = (15) 
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FIG. 2: Structure of axisymmetric equilibrium von Karman flow driven by electro- magnetic system 
at the boundary for Mach number M = 1, fluid Reynolds number Re = 200 and magnetic Reynolds 
number Rm = 20: (a) number density; (b) velocity; (c) magnetic field. Cross-sections in r — z plane 
are given. Left panels of (b) and (c) show stream lines of poloidal parts (r- and z-components) of 
flux nv and magnetic field b, respectively, superimposed on absolute values of these parts depicted 
in colors. Right panels of (b) and (c) show azirQuthal components of corresponding fields. 



(perfectly conducting walls), and the velocity conditions are: 

v| r= i = ze v , (16) 
v|*=_i = -re v , 
v| 2= i = re v 

(no-slip differentially rotating rigid walls). These velocity boundary conditions correspond 
approximately to von Karman flow: the top and bottom end-caps are counter-rotating with 
the same angular velocity, and the side wall has a linear dependence of angular velocity on 
z to match the rotation of top and bottom end-caps. This flow is the primary object of our 
dynamo study. 

We briefly remark on the spatial and temporal resolution used in these simulations. 
For spatial discretization NIMROD employs a high order finite element method in r- and 
z-directions and a pseudospectral method in periodic (^-direction with dependence e tmip 
for each Fourier harmonic (integer m represents the azimuthal mode number). The basis 
functions of the finite elements are polynomials. In all of the simulations presented here we 
have used a uniform meshing of the r — z plane with 8x16 finite elements each of polynomial 
degree 3, and 11 Fourier harmonics in the azimuthal (^-direction. This spatial resolution 
appears to be sufficient for the laminar flows under consideration. To verify the simulation 
results obtained at this resolution, we have repeated a nonlinear MHD run introduced in 
section V using mesh with 16 x 32 finite elements in r — z plane and 11 azimuthal modes. 
The time dynamics of the flow and the magnetic field were fully reproduced, which confirms 
that a converged solution is obtained already at the coarser resolution. The solutions are 
time-evolved using a semi-implicit staggered leap-frogging algorithm, which is fully detailed 
in Refs. [16, 17]. We should emphasize here that the algorithm employed in NIMROD 
can be made numerically stable for arbitrarily large time-steps in both single-fluid and Hall 
MHD by choosing the appropriate coefficients of the semi-implicit operators [17]. However, 
in order to accurately model the temporal behavior of the system with significant flows, in 
the present simulations we have used an adaptive time-step based on the explicit Courant- 
Friedrichs-Lewy (CFL) condition for advection. Even in the dynamo simulations, we find 
that advection by the axisymmetric velocity field dominates the CFL criteria, rather than 
either the Alfven or whistler waves associated with the relatively weak magnetic fields. 
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III. HYDRO DYNAMICAL EQUILIBRIUM AND STABILITY 



In this section we consider the basic hydrodynamical (no magnetic field, b = 0) properties 
of von Karman flow in cylinder, in which the plasma is stirred at the edge via the boundary 
conditions of set II (16). We start with axisymmetric case, assuming that physical quantities 
do not depend on ip, i.e., d/dip = 0. Axisymmetric equilibrium flow structure for fluid 
Reynolds number Re = 200 and Mach number M — 1 is shown in Fig. 3. The velocity 
components are either symmetric (radial v r ) or antisymmetric (azimuthal v v and axial v z ) 
with respect to equatorial plane z = 0, and the azimuthal velocity has maxima at the corners 
of the cylinder. The flow develops two cells of poloidal (in the r — z plane) circulation with 
inward direction at the equatorial plane. Such flow pattern leads to the stratification of 
density, which builds up near the corners of the cylinder. 

The axisymmetric equilibrium von Karman flow becomes unstable with respect to non- 
axisymmetric perturbations when the fluid Reynolds exceeds a critical value (Fig. 4). This 
instability is the Kelvin-Helmholtz type instability, occurring due to the presence of unstable 
velocity shear in the fluid. For Mach number M = 1 the critical value of Reynolds number 
is about Re ~ 160. As shown in Fig. 4, within approximately one viscous time the unstable 
modes grow and saturate at a new equilibrium state, which consists of the axisymmetric 
part and relatively small non-axisymmetric distortions with even azimuthal mode numbers, 
m = 2, 4, 6, . . .. Such transition to non-axisymmetric equilibrium is crucial for the dynamos 
considered here. 

The structure of the equilibrium von Karman flow, in particular the amplitude of the non- 
axisymmetric distortions, depends on Reynolds and Mach numbers (Fig. 5). As we discuss 
in section IV, such dependence affects the critical value of magnetic Reynolds number above 
which the dynamo is excited. 

IV. KINEMATIC DYNAMO 

Our first step in this dynamo study is to solve the kinematic dynamo problem: determin- 
ing the possibility of self-generation of magnetic field for a given flow structure. For fixed 
fluid Reynolds and Mach numbers we solve the stationary time independent (d/dr = 0) 
continuity (10) and Navier-Stokes equations (11) and find the steady-state hydrodynamic 
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FIG. 3: Structure of axisymmetric equilibrium von Karman flow driven by differentially rotating 
walls for Mach number M = 1 and fluid Reynolds number Re = 200: (a) number density; (b) 
velocity. The same elements as in Fig. 2 are shown. 

equilibrium velocity v cq (which includes possible non-axisymmetric distortions). Using that 
velocity, we solve induction equation (12) as an eigenvalue problem for magnetic field with 
different magnetic Reynolds numbers 

1 



7 b = V x (v cq x b) + ^V 2 b, 



(17) 



where 7 is an eigenvalue. This allows us to determine the critical magnetic Reynolds above 
which the dynamo excitation is possible. The results are presented in Fig. 6. The dependence 
of the critical magnetic Reynolds on the fluid Reynolds for Mach number M — 1 is shown 
in Fig. 6(a). The vertical line at Re ~ 160 separates the regions of axisymmetric (Re < 
160) and non-axisymmetric (Re > 160) von Karman flows. Our simulations show that 
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FIG. 4: Time dynamics of kinetic energy of different azimuthal harmonics in purely hydrodynamical 
(b = 0) von Karman flow for Mach number M = 1 and fluid Reynolds number Re = 200. 
Corresponding azimuthal mode numbers m are shown. 

the kinematic dynamos are not possible in axisymmetric flows. In non-axisymmetric flow at 
sufficiently high magnetic Reynolds number Rm, the dynamo appears as a growing magnetic 
field with odd azimuthal harmonics, m = 1, 3, 5, . . .. 

In our study of the plasma dynamo we use a compressible fluid model. Compressibility is 
related to the Mach number (4) - the ratio of the peak driving velocity to the sound speed. 
In general, the higher Mach number the more compressible the fluid, i.e., the more stratified 
its density. An increase in the Mach number leads to changes in the equilibrium velocity 
field v eq and, in particular, decreases the energy in the non-axisymmetric components of the 
flow (Fig. 5(b)). This, in turn, affects the kinematic dynamo problem (17) by increasing the 
value of the critical Rm (Fig. 6(b)). 

From a comparison of Figs. 5 and 6 we can conclude that the dynamo action in the 
case under consideration is related to the non-axisymmetric part of the von Karman flow: 
the stronger the non-axisymmetric distortions, the lower the threshold value of magnetic 
Reynolds number for the onset of the dynamo. 

V. NON-LINEAR SATURATION OF DYNAMO FIELD 

Next we analyze the state of the fully saturated dynamo generated magnetic field and its 
back reaction on the flow. In this section we report the results of fully non-linear simulations 
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FIG. 5: Kinetic energy of different azimuthal harmonics in hydrodynamically stable von Karman 
flow as a function of (a) fluid Reynolds number Re for Mach number M = 1; (b) Mach number M 
for fluid Reynolds number Re = 200. Corresponding azimuthal mode numbers m are shown. 

of system (10)-(12). The non-dimensional parameters are chosen to be Re = 200, Rm = 
1000, M = 1 and e = 0.05 — 1.0 (for Hall MHD runs); simulations with these parameters 
achieve fully saturated laminar states for the dynamo fields and the fluid flows. 

Figs. 7(a,b) demonstrate the time dynamics of the kinetic and magnetic energies for a 
single-fluid MHD case (e = 0). After the initial transient phase (r w 100), the flow becomes 
stationary. It is primarily axisymmetric, with non-axisymmetric distortion consisting of 
even azimuthal harmonics (m = 2, 4, 6,...) and containing only about 1% of the total 
kinetic energy. In such a flow the dynamo is excited: magnetic field of odd harmonics 
(m = 1, 3, 5, . . .) grows exponentially in time until it saturates at r rj 1500. In some 
sense, the saturated magnetic field is in equipartition with the non-axisymmetric part of 
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FIG. 6: Critical magnetic Reynolds number Rm as a function of (a) fluid Reynolds number Re for 
Mach number M = 1; (b) Mach number M for fluid Reynolds number Re = 200. Vertical line in 
(a) separates the regions of axisymmetric and non-axisymmetric equilibrium von Karman flows. 

the flow. Here -E mag is about 0.3% of the total energy of the flow -E^m- Due to the lack of 
an axisymmetric magnetic field (m = 0) and the small amplitude of the non-axisymmetric 
fields, the back reaction of the dynamo magnetic field on the axisymmetric flow is very weak 
and the imposed von Karman flow is essentially unmodified. The structure of the saturated 
magnetic field, with the m = 1 azimuthal harmonic visibly dominating the overall structure, 
is shown in Fig. 8(a). 

In the Hall MHD case (e = 0.05) the dynamics of the system is different (Figs. 7(c,d)). 
After the non-axisymmetric flow with even modes develops (at r w 100) the odd harmonics 
of the magnetic field start to grow, as in single-fluid MHD case. Now however, when the odd 
harmonics of the magnetic field become large enough (at r pa 500), the Hall effect becomes 
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FIG. 7: Time dynamics of kinetic E^ n and magnetic E m£Lg energies of different azimuthal modes 
in von Karman flow for Mach number M = 1, fluid Reynolds Re = 200, and magnetic Reynolds 
Rm = 1000, with azimuthal mode numbers m labelled. (a,b) Single-fluid MHD case (e = 0). Flows 
are of even m modes while fields are odd in m. (c,d) Hall MHD case (e = 0.05). Initial behavior 
is similar to the single-fluid MHD case, but as the fields become strong (r ~ 500), the Hall effect 
becomes important. The final equilibrium includes both odd and even m. 



significant and other harmonics of the flow and the magnetic field grow, breaking the initial 
even-odd symmetry. The saturated state of the flow and the magnetic field contains all 
harmonics, though the dominating parts are similar to MHD case. The structure of the 
saturated Hall dynamo field for e = 0.5 is shown in Fig. 8(b). 

Fig. 9 shows the dependence of the magnetic energy for different azimuthal modes in 
the saturated state on the Hall number e. It is worth noting that the presence of the Hall 
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(a) Single-fluid MHD dynamo, 8=0 
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(b) Hall MHD dynamo, 8=0.5 




FIG. 8: Magnetic field lines of saturated dynamos: (a) single-fluid MHD case (e = 0); (b) Hall 
MHD case (e = 0.5). Thickness of the line is proportional to the magnitude of the field, while red 
(blue) color corresponds to upward (downward) direction of the field. 
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effect not only breaks the even-odd symmetry of the system, but also reduces the energy of 
the saturated dynamo field. For e > 0.2 the energy of the saturated dynamo field scales as 
-E'mag oc e~ 2 (b oc e^ 1 for the field amplitude). This is a direct consequence of the magnetic 
induction equation (12), which in a saturated state (d/dr = 0) now reads: 

Vx (vxb--(Vxb) xb) +— V 2 b = 0. (18) 
\ n J Rm 

Indeed, if the magnetic Reynolds number Rm and functions v and n are independent of e, 
then the solution for the magnetic field b can be written as 

b = |, (19) 

where h is a vector-function independent of e satisfying equation 

V x ( v x h - -(V x h) x h ) + — V 2 h = (20) 
\ n J Rm 

and boundary conditions (15). In the case under consideration it is clear that the scaling used 
in (19) is valid only asymptotically for large Hall numbers e when the magnetic field b is small 
and its influence on the flow is negligible. Under these assumptions, the plasma velocity v 
and density n do not depend on e and correspond to the purely hydrodynamical equilibrium 
that is determined by the fluid Reynolds Re and Mach M numbers (as described in section 
III). Thus, for relatively large Hall parameters (e > 0.2) the saturation of the dynamo takes 
place due to the back reaction of the Hall term in induction equation (12) long before the 
amplitude of the magnetic field is large enough to change the flow significantly. This is in 
contrast to the single-fluid MHD dynamo (and the Hall dynamo with e << 0.1), where the 
saturated state is achieved due to modification of the flow profile by the growing dynamo 
field before the Hall effect plays a considerable role. Scaling (19) suggests that the Hall 
effect is unfavorable for the dynamo. 

We note Refs. [22-24], where results of simulations for a single-fluid MHD dynamo in 
a sphere with counter-rotating incompressible flows are reported. Despite the difference in 
geometry, our MHD results are in qualitative agreement with the cases of laminar dynamos 
from Refs. [22, 23], which have dominant m = and m = 2 harmonics in kinetic energy and 
m = 1 harmonic in magnetic energy during the saturation phase. However, in contrast to 
results of Ref. [24], we do not observe the generation of an axial magnetic dipole with m = 
and quasi-periodic oscillations. 
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FIG. 9: Magnetic energy of different azimuthal harmonics in saturated Hall MHD dynamo as a 
function of the Hall number e. Azimuthal mode numbers m are shown. Dashed line corresponds 
to scaling -E mag oc e -2 . 

VI. DISCUSSION 

We have performed numerical simulations of a plasma dynamo for cylindrical von Karman 
flow using the extended MHD code NIMROD. These NIMROD simulations provide numer- 
ical support for the Madison Plasma Couette Experiment (MPCX). We have demonstrated 
that sustained dynamo action and self-generation of magnetic field can be attained for pa- 
rameters that are achievable in the experiment. Our results show that the critical magnetic 
Reynolds number required for dynamo excitation strongly depends on the plasma compress- 
ibility: the more compressible the fluid, the higher the critical Rm (Fig. 6). Inclusion of 
two-fluid effects into the model (in the form of the Hall term in the induction equation) does 
not influence the critical Rm, but does change the structure of the saturated flow and the 
dynamo field (Fig. 7). The effect of the Hall term on dynamo field is negative: the energy 
of the saturated magnetic field scales as 1/e 2 when the Hall number is e > 0.2 (Fig. 9). 

The simulations show that the presence of non-axisymmetric distortions in the flow plays 
a decisive role in dynamo excitation. Such distortions with even azimuthal mode numbers 
(m = 2, 4, 6, . . .) appear in the flow only for modest values of fluid Reynolds number 
Re > 160 when the axisymmetric shear flow becomes hydrodynamically unstable. Therefore, 
in order to observe the dynamo effect in the MPCX, the plasma has to be driven above a 
critical hydrodynamical threshold. For a helium plasma with number density n = 10 18 
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m~ 3 , electron temperature T e = 16 eV and ion temperature Tj = 0.9 eV, the peak driving 
velocity should be V = 20 km/s. 

Another issue is related to the detection of the laminar dynamo field in the experiment. 
Since the saturated magnetic energy of the dynamo field is only a small fraction of the 
total kinetic energy, the resulting dynamo field is relatively weak. For the helium plasma 
considered above, the volume-averaged saturated dynamo magnetic field is B Q ps 0.1 Gauss. 
Such a small field is still detectable even on the background of the much stronger multi- 
cusp field from the rings of magnets (about 10 3 Gauss near the walls) due to the different 
azimuthal symmetries of these two fields. 

Lastly, we remark on the model used in our simulations. The single-fluid MHD model 
(e = 0) is adequate for predicting critical magnetic Reynolds numbers and the thresholds 
for sustained dynamo action. The fully non-linear saturated dynamo state differs however 
when the Hall effect is included. It appears that the isothermal Hall MHD model (e > 0) is a 
good "rough" approximation for the plasma under experimental conditions, but this model 
also does not capture the full details of the plasma dynamics. Other effects such as thermal 
conductivity, electron pressure, and anisotropic viscosity can all play important roles in a 
real plasma experiment. These effects will be the subject of future studies. 
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